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We derive expressions for various correlators of the gauge field and find the propagators in 
Hamiltonian dynamics which employs a new gauge — 0. This gauge is a part of the wedge 
form of relativistic dynamics suggested earlier as a tool for the study of quantum dynamics in ultra- 
relativistic heavy ion collisions. We prove that the gauge is completely fixed. The gauge field is 
quantized and the field of radiation and the longitudinal fields are unambiguously separated. The 
new gauge puts the quark and gluon fields of the colliding hadrons in one Hilbert space and thus 
allows one to avoid factorization. 
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In two previous papers of this cycle further quoted as papers [I] and [II]), we explained the physical motivation 
of the "wedge form of dynamics" as a promising tool to explore the processes which take place during the collision 
of two heavy ions. In compliance with the general definition of dynamics given by Dirac 1^, the wedge form includes 
, \ a specific definition of the quantum mechanical observables on space-like surfaces, as well as the means to describe 
I> [ the evolution of the observables from an "earlier" space-like surface to a "later" one. Unlike the other forms, the 
I wedge form explicitly refers to the two main geometrical features of the phenomenon, i.e., the strong localization of 
I the initial interaction and, as a consequence, the absence of translational invariance in the temporal and longitudinal 
' directions. 



In the wedge form of dynamics, the states of the quark and gluon fields are defined on the space-like hyper-surfaces 
of the constant proper time r, — — z'^. The states of fermion fields were discussed in paper [II]. In this paper, we 
Oh continue the study of the gluons and augment our previous consideration by the gauge condition — 0. This simple 
O i' idea solves several problems. First, this gauge condition is boost-invariant and thus complies with the symmetry of 
pH ' the collision. Second, it becomes possible to treat two different light-front gauges (which describe gluons from each 
nucleus of the initial state separately) as the two limits of this single gauge. Therefore, the new approach keeps 
important connections with the theory of deep inelastic ep-scattering (DIS). This fact is vital for the subsequent 
^ ' calculations since e — p DIS is the only existing source of data on nucleon structure in high-energy collisions. The 
■ ■ ■ ' approach based on quantum field kinetics (QFK) allows one to treat both the nuclear collision and e — p DIS as the 
similar transient processes. Third, after the collision, this kind of gauge becomes a local temporal axial gauge, thus 
providing a smooth transition to the Bjorken regime of the boost-invariant expansion. 

Most of this paper is technical, and any relevant physical discussion of the results appears only after their mathe- 
matical derivation. These results were summarized in paper [II]. Since the first interaction of two finite-sized nuclei 
is strongly localized, the geometrical symmetry of the final state is manifestly broken and the observables of wedge 
dynamics are essentially defined on the curved space-like hypersurface. For the fermion field this have led to an 
obvious Thomas precession. Similar orientation effects happen in the case of the vector field also. The dynamics of 
the gauge field is rich, and the procedure of its quantization triggers many puzzles that can be traced back to the 
classical roots of the gauge field theory. 



In section II A, we derive equations of motion for the gauge field in the gauge — 0, find the Hamiltonian variables 



and the normalization condition. The equations of motion are linearized and the modes of the free radiation field are 
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obtained in section 11 B. In section p C| , the retarded propagator of the perturbation theory is found as the response 
function of the field on the external current. This part of the calculation turned out to be the most durable, since the 
gauge condition is inhomogeneous and none of the modern methods is effective. However, the old-fashioned variation 
of parameters does work. The most important result of this paper, separation the transverse and longitudinal parts 
of the gluon propagator is obtained here. We essentially base calculation of the quark self-energy in expanding quark- 



gluon system on this result. These calculations are presented in the next paper. In section II D, we show that the 
previously obtained propagator solves the initial data problem for the gauge field. Unlike in the homogeneous axial 
gauges, the propagators of the gauge — do not have any spurious poles. 



Section [II is devoted to the quantization of the vector field in the gauge — 0. We begin in III A with the 
proof of the fact that the gauge — can be completely fixed provided the physical charge density, rj^ vanishes at 
T — 0, the moment of the first touch of the nuclei. This is exactly what can be expected from the nuclei color lessness. 
Then, the Gauss law can be unquestionably used to eliminate the unphysical degrees of freedom in the equations of 



motion. We continue in III B with a computation of the Wightman functions, and study the causal properties of the 



commutators in |III C| . The latter appears to be abnormal; the Riemann function is not symmetric and penetrates the 
exterior of the light cone. However, the behavior of the observables is fully causal and the procedure of the canonical 



quantization is accomplished in HI D . Even though it is impossible to introduce transverse and longitudinal currents 
(as it is customary for the homogeneous gauge conditions) and thus fully separate the dynamics of the corresponding 
fields, we found it useful to discriminate the various field patterns by the type of their propagation. The propagator 
of the transverse field is sensitive to the light cone boundaries while the longitudinal and instantaneous parts of the 



field do not propagate. These coordinate form of these two fragments of the response function is derived in section IV. 
Ultimately, the longitudinal part of the gluon propagator appeared to be of the greatest importance for the dynamics 
of the screening effects at the early stage of the collision. 

In Appendix 4, we study the limiting behavior of the propagator in the central rapidity region and in the vicinity 
of the null-planes and show that propagators of the gauges — and = are recovered. It is important 
that the spurious poles are recovered only in the unphysical limit of infinite rapidity. This result is of practical 
importance because it establishes the connection between the new approach and the existing theory of the deeply 
inelastic processes at high energies. 



II. THE CLASSICAL TREATMENT 



The final goal of this paper is to build a quantum theory of the vector gauge field in the expanding geometry of 
nuclear collision. Development of a quantum theory always begins with its classical counterpart which provides the 
one-particle wave functions (which later serve as quantum states) and the classical Green functions (which later become 
the propagators of quantum theory). Furthermore, the quantum propagator of gauge field includes the longitudinal 
part which can be found only by classical analysis. The classical part of this program is the subject of this section. 



A. Classical equations of motion 



We now consider the case of pure glue-dynamics. We denote the gluon field in the fundamental representation of 
the color group as ^^(x) — t°'A'^{x). Gonsequently, we have the field tensor, 

i^M- = = ^^M^^ - 25-^^ = df^A, - d,A^ - ig[Af,, A,], 
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where = 9^ — ig[Af^{x), ...] is the covariant derivative on the local color group. The gauge invariant action of the 
theory looks as follows, 

S = J C{x)d^x = j [-\g^^\x)g''%x)F^,{x)F^,{x) - fA^]^d^x. (2.1) 

Its variation with respect to the gluon field yields the Lagrangian equations of motion, 

dx[{~gf'h''h''"F,,] - ^g{-gf'^[A^,g^Y"F,^] = {-gf/^f , (2.2) 

where is the color current of the fermion fields and g = det|g^^|. The equations are twice covariant, i.e., with 
respect to the gauge transformations in color space and the arbitrary transformations of the coordinates. In what 
follows, we shall employ the special coordinates associated with the constant proper time hyper-surfaces inside the 
light cone of the collision point t = z = 0. The new coordinates parameterize the Minkowski coordinates (t, x, y, z) as 
(t coshry, X, ?/, r sinh ?/). In addition, we impose the gauge condition Ar = 0. The corresponding gauge transformation 
is well defined. Indeed, let A^(x) be an arbitrary field configuration and A'^{x) its gauge transform with the generator 



J7(r, 77,rl) = P^exp{- / A^(r', 77, rDdr'}, 
Jo 



(2.3) 



then the new field, A'^ — U A^U ^ + d^UU ^, obeys the condition — 0. Imposing this gauge condition we arrive 
at the system of four equations: 

C{x) = -dndrA„ + TdrdrAr - ig\^\A„,drA„\ + T\Ar,drA,\} - rf = 0, (2.4) 
T r 

^ drTdrAr + —d^{d,^Ar — + Tds{dsAr — drAg) 

-ig{-d,j[A^,Ar] + Tds[A,,Ar] + -[^^, F,^r] + T[As,F,r]} - rf = , (2.5) 

T r 

- dr-drA,, + -dr{drA,, - d^Ar) - ig[-dr[Ar , A,,] + -[Ar, F„,]] - t/' = (2.6) 
T T T T 

Here, we use the Latin indices from r to w for the transverse x- and y-components (r, w = 1,2). We shall also use 
the arrows over the letters to denote the two-dimensional vectors, e.g., k = (k^, ky), \k\ = k±. The Latin indices from 
i to n {i, n = 1,2, 3) will be used for the three-dimensional internal coordinates = (x, y, j]) on the hyper-surface 
T = const. The metric tensor has only diagonal components grr — ^gxx ~ ^Syy — 1j Svn ~ ■ first of these 



equations (2.4) contains no second order time derivatives and is a constraint rather than a dynamical equation. The 
constraint weakly equals to zero in classical Hamiltonian dynamics and serves as a condition imposed on physical 
states in the quantum theory. The canonical momenta of the theory are as follows: 

^-Fr^ = -A,j, Pr=rFrr=TAr ■ (2.7) 
T T 

Hereafter, the dot above the letter denotes a derivative with respect to the Hamiltonian time r. Because of the 
gauge condition, the canonical momenta do not contain the color commutators. After excluding the velocities, the 
Hamiltonian can be written down in the canonical variables, 

H = J drjdrlT{^py'^ + ^fp' + ^F^rF^r + ^Fr^Frs + fA„ + fAr} . (2.8) 
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Then the equations (2.5) and ( |2.6| ) are immediately recognized as the Hamiltonian equations of motion. The Poisson 
bracket of the constraint C with the Hamiltonian vanishes, thus creating the generator of the residual gauge trans- 
formations which are tangent to the hyper-surface. Conservation of the constraint is a direct consequence of the 
Lagrange (or Hamiltonian) classical equations of motion as well. 

The normalization condition for the one-particle solutions is obviously derived from the charge conservation law. 
For the gauge field, this is impossible. Therefore, we shall accept the condition which supports self-adjointness of the 
homogeneous system after its linearization. This leads to a natural definition for the scalar product of the states of 
the vector field in the gauge — , 

{V, W)= I df] f d^rTg'''V*idrWk , (2.9) 



where g**"' is the metric tensor of the three-dimensional internal geometry of the hyper-surface r = const . This norm 
of the one-particle states prevents them from flowing out of the interior of the past and future light wedges of the 
interaction plane. 

B. Modes of the free radiation field. Field of the static source 

As a tool for the future development of the perturbation theory, we need to find the propagators and Wightman 
functions when the nonlinear self-interaction of the gluon field is switched off. In this case, the system of equations 
for the non- vanishing components of the vector potential and the constraint look as follows, 

[drTdr - -dl - ra2]A, + dr\TdsA, + -d,,A.^\ - -rf , (2.10) 

T T 

[0^}Ldr--d',]A, + -d,d,A, = -rf , (2.11) 

C{x) = -d„drA„ + TdrdrAr - rf = , (2.12) 
r 

where j'^ includes all kinds of the color currents. An explicit form of the solution for the homogeneous system is found 
in Appendix 1. In compliance with the gauge condition (which explicitly eliminates one of four field components) we 
find three modes V'-'^'' of the free vector field. Two transverse modes obey Gauss' law without the charge and have 



the unit norm (see Appendix 1) with respect to the scalar product (2.9) 



vkyRl^X-^Ak^T) leW^'v (2.13) 



V -4'^,.(fc±r) 



The mode F^^^ is constructed from the functions R^J^_i^{k_\_T) — i?^^Ljj,(fcj_r|s) corresponding to the boundary 
condition of vanishing gauge field at r = 0. This guarantees continuous behavior of the field at r = 0. Indeed, as 
T ^ 0, the normal and the tangential directions become degenerate. As long a.s A'^ — is the gauge condition, 
continuity requires that A'' as r ^ 0. 

It is instructive to know the physical components of the electric and magnetic fields of these modes, — 
V=^g"" An and 6™ ^ -(2V^)-ie"'"Fi„: 
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f ('^"(x) = = ( l^t 1 i?^'Ufc±^)e"''+'^'' ■ (2.14) 

The mode F^^) 

can be obtained from the mode V^^^ by a simple interchange of its electric and magnetic fields. Using 
standard wave-guide terminology, one may call mode y*^^^ as the "transverse electric mode" and the mode V'^'^^ as the 



"transverse magnetic mode". Equations ( 2.14 ) indicate, that the field strength tensor of the free radiation field obeys 
the condition, (F*)* = —F. Therefore, certain linear combinations of the modes V'^^^ and V'^'^^ may be analytically 
continued to Euclidean space as self-dual solutions of the field equations. 

An equivalent full set of the transverse modes carries (instead of the boost v) the quantum number 6, (rapidity), 
i.e. ko — k± cosh^?, = k±_ sinh^?. These functions can be obtained by means of the Fourier transform. 



and have the following form. 




where 



= [ ) e--— .g(x) = I kj, 1 , (2.16) 



/i(t,?7) = k± smh{e - 77) / e-'''^^' "'"'^^'^-"Ut' = itanh{9 - ^)(e-'fc-L^™>'h(e-'?) _ l) 







^—ik±T cosh{0~ri) -t ~ik±T cosh{9 — ri) 

h{r,v) = kl e-^'^.rcosHO-,)^,^^^ —+^k^T — ^ . (2.17) 

cosh"' (6* - 77) cosh(6'-?7) 







The norm of the Coulomb mode V'^'^\ as defined by Eq. (2.£), equals zero, and it is orthogonal to F^^^^ and V'^ 



(2) 



Though this solution obeys the equations of motion without the current, it does not obey Gauss' law without a charge. 
Therefore, it should be discarded in the decomposition of the radiation field. However, it should be kept if we consider 
the radiation field in the presence of a static source with the r-independent density p(fc, v) — rjZ (r) — const{T). In 
this case, its definition can be completed using Gauss' law: 

k.i'^ ' {2iTYik\ \ V Qi^i„[k±_T) J ^ ' 

The coordinate form of this solution is noteworthy. The physical components, = gg"''A; , of the electric field 
of the "r-static" source can be written in the integral form. 



^lstat)i'^^^i^Vi) = j dr2dinKi{T]ri - f2,m - m)p{r2,m) , (2.19) 

with the kernel 

K(r-rn)- f^^—^f >'rS^Mk±r) \ Ojr - r^) f r coshr^jd/dxn \ 1 

lU(T,r,fj)-j ^^^^3 [i^kls^iMk^r) J - An \ 9/a(rsinhry) )r^^ ' ^^'^^^ 

where R12 = {r\ + sinh^ riY^'^ = [{ri — ^2)^ + sinh^(r;i — 772)]"'^^^ , is the distance between the points (ri, 771) 
and (7^2, '72) iu the internal geometry of the surface r = const. The technical details of the derivation of the last 
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expression will be presented in Sec. IV. Eq. (2.20) is an analog of Coulomb's law of electrostatics, except that now the 
source has a density which is static with respect to the Hamiltonian time r. In fact, the source is static if it expands 
in such a way that its physical component J"^ — rj'^ [t, rj, r) does not depend on t. These expressions will be helpful 
in recognizing the origin of various terms in the full propagator which is calculated below. 

C. Propagator in the gauge = 

The calculation of the propagator in the gauge — (associated with the system of the curved surfaces r — const) 
meets several problems. Three methods are commonly used in field theory. One of them strongly appeals to the Fourier 
analysis in the plane Minkowski space which is not applicable now because the metric itself is coordinate-dependent. 
The second method uses the path-integral formulation which is also ineffective because of the explicit coordinate 
dependence of the gauge-fixing term in the Lagrangian. One could also try to study the spectrum of the matrix 
differential operator, to find its eigen-functions, and to use the standard expression for the resolvent. However, the 
extension of the system for the non-zero eigenvalues leads to unwieldy equations. On the other hand, the Green 
function of the perturbation theory must coincide with the one which solves the problem of the gauge field interaction 
with the classical "external" current. For this reason, we shall compute the Green function in a most straightforward 
way; we shall look for the partial solution of the inhomogeneous system using the old-fashioned method of "variation 
of parameters" . This method will immediately separate the radiation and the longitudinal parts of the retarded 
propagator. All other methods would require an additional analysis for this purpose. 

Let us start the derivation of the propagator in the gauge = by obtaining the separate differential equations 

for the ry-component of the magnetic field, "if = dyA^ — dxAy, the transverse divergence of the electric field, tp — 

• • • 

T{dxAx + ^yAy)^ and the 77-component of the electric field, a = Ati/t ■ In terms of the Fourier components with 

respect to the spatial coordinates, these equations read as 

[^r + ^9. + ^+fci]%.M=-j'^(fc,^,T) , (2.21) 

[drTdr + —](f{k,v,T)^iTvkia{k,v,T)^-dr[T'^f(k,i',T)] , (2.22) 
r 



[drTdr + Tkl]a{k,iy.T) (^(fc, z^, t) = -9, [t^j" (fc, I/, r)] , (2.23) 

T 

where j"^ — dyjx ~ dxjy, — dxjx + dyjy. Using the constraint conservation, which may be explicitly integrated to 
tp{k, t)) + iiya{k,v,T) — Tj'^ {k,iy,T) — —pQ{k,u) — coiistir) , (2.24) 
one easily obtains two independent equations for ip{k, v, r) and a(fc, v, r): 

[dl + -dr+'^ + klMk, V, t) = fci \p{k, V, t) - v)] - -dAT^fik, V, r)) = r , (2.25) 

T T 

[9^ + is, + 4 + kl]a{k, V, t) = ^[p(fc, V, r) - po(fc, v)\ - -dr{T^f{k, v, r)) = /" , (2.26) 

The constant of integration po(^i ha-s the meaning of the arbitrary static charge density and it should be retained 
until Gauss' law is explicitly imposed on the solution. In what follows, we shall not write it explicitly, keeping in mind 
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that it is included in the true charge density p(k, v, t). Since wedge dynamics has a selected time moment r = 0, the 
constant of integration po{k, v) can be associated with the initial data, namely, with the charge density at r = 0. As 
we shall see soon, a proper choice of po will be needed in order to fix the gauge = Q completely. 
Equations ( ^.2lD , ( |2.25| ) and ( p^ ) can be solved by the method of "variation of parameters" : 
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TT? 

^(r)--/ T2dT2n{T,T2)f{T2) , (2.27) 







where T stands for any one of the unknown functions in these equations, and / for the corresponding right hand side. 
The kernel 

is the usual bilinear form which is built from the linearly-independent solutions of the homogeneous equation. [ The 
Wronskian of these solutions is exactly 4/i7rr2. ] Taking J- = 'f, we obtain the first equation for the components 
Ax{k, v, t) and Ay{k, v, r) of the vector potential: 

^{k,V,Tl) = l[-kyA^+krAy] = — T2dT-2'H{Tx,T2) i[-kyf{T2) + k^jy{T2)] . (2.28) 

4 Jo 

In order to find the second equation for the x- and y-components and the equation for A,-i{k, v, r), we must integrate 
twice, i.e., 

in dr' r' ^ 1 ^ 

<^{k,V,Tl) = l[k^A,+kyAy\^— — ■H{T',T2)T2dT2[~kl p{k,V,T2) + —dr^{T^f{k,V,T2))] , (2.29) 

4 Jo '''Jo '^2 

77r /""^^ C^' iv -> 1 

A^{k, V, Tl) = - / r'dr' / H(t', T2)T2dT2 {— p{k, V, T2) + -9., (r|j''(fc, V, T2))] . (2.30) 
4 Jo Jo ^2 "^2 

The integration over T2 recovers the electric fields at the moment r', whilst the integration over r' gives the vector 
potential at the moment ri. It is convenient to start with the second of these integrations which has the limits 
T2 < t' < Tl. Let us consider the main line of the calculations in detail, using the 77-component as an example. The 



first integration follows the formula (Bl), 



fc'^+i r{TrH\^){kiT')d^' = i?i':L(fc±n) - Bli„,{k^r2) , (2.31) 

J T2 

and the terms emerging from the lower limit T2 can be conveniently transformed according to the relation (see 
Appendix 2), 

R%{k^r2)H\l\k^T2) - Rf,,Akl_T2)H\l\k^T2) = ^S^,,Ak^T2) . (2.32) 
As a result, one obtains, e.g., the following formula for A,j(fc, v, ti): 

Mk,l^,n) = r T2dT2[R['iik^T,)Hl^\k^T2) ~ ^1 (fc^n )i? (fc^T^ ) - 1 Si,,;, (fc^Ts )] 

x[^p{k,l^,T2) + -dr,{Tif{k,V,T2))] . (2.33) 

In order to eliminate the charge density p from the integrand and to separate the transverse and the longitudinal parts 
of the propagator, all the terms of this formula should be integrated by parts, explicitly accounting for the charge 
conservation, which reads as 
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iT[k^f{k,iy,T) + k^f{k,iy,T) + uf{k,v,T)] +drp{k,v,T) = 



(2.34) 



We have in sequence: 



Jo 



^'^pi^r,)H[^{k^r,) = ^v f 

Jo "^2 Jo 



p{T2)dT2 



(2.35) 



IV dT2H^,p{k^T,)drArb'{r2))=r!f{T,)H^,l\k^r,)+ / T,dT2[R[%{k^T,) + R^ll^ik^r^)]^ {r^). (2.36) 
In a similar way we have, 



dQ-l,iu{k^T2) 



p{T2)dT2 



lU I p{T2)si_^i,{k±_T2) ^ iv , , 

lo T2 Jo dT2 

il'Q-l,iu{k±Ti)p{Ti) - V ( T2dT2Q-l^ru{k±T2)[k^f{T2) + kyf{T2) + l'j''(T2)] 
^0 



(2.37) 



dT2S,,Uk±T2)dr,{Tif{T2))^Tif{n)siMk±n) + i^^ / Tadra [Q- 1,„ (fc^ ) - Qi,., (fc^Ta )b-'' (t2 ) . (2.38) 
/o Jo 

Assembhng these pieces together and repeating the same calculations for the function $ one obtains three different 

terms which contribute to the field A produced by the current j; A = A*^*""^ + A'^^ + ^(*"^*). 



The transverse field A^*''^ is defined by the integral terms in the R.H.S. of Eqs. (2.35) and (2.36). It can be 
conveniently written down in the following form, 



Af''^(xi)= / d^X2e{Ti~T2)A[Z\xi,X2)r{x2) 



where 



A 



r\x,y) ^-i dv d^k y [V^y (x)V^y* (y) - V^y*(x)V'^i^ (y)] 



(2.39) 



(2.40) 



A=l,2 



which can be easily recognized as the Riemann function of the original homogeneous hyperbolic system. The Riemann 
function solves the boundary value problem for the evolution of the free radiation field. It is obtained immediately 



in the form of the bilinear expansion over the full set of solutions (2.13) of the homogeneous system. In fact, this is 
a sole evidence that A'^*''^ may be associated with the transverse part of the propagator. Then the remaining part is 
the propagator (response function) for the longitudinal field. 



The dynamical longitudinal field A^^' originates from the integral terms in the R.H.S. of Eqs. (2.37) and (2.38): 
A|^^(Ti,r;i,fi) = ^ ' T2dT2 J dTi2d^f2A\^{T2;m " m,ri - r2)r{r2,V2,r2) ■ (2.41) 
The kernel of this representation. 



aJ^'(t2;7?i -r?2,ri -fi) 



r diyd^k 






1 i^n^kl 


V 


1 



ksQ-l,iu{k±T2) 
vQl,ii,{k^T2) 



(2.42) 



does not allow for the bilinear expansion with two temporal arguments, and, as we shall see in a while, the retarded 



character of the integration in Eq. (2.41) is not sensitive to the light cone boundaries. In fact, the electric field 
El — Ai is simultaneous with the current . 
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The instantaneous part of the solution conies from the boundary terms in Eqs. ( ^.35| )-(2.38) which were generated 
via integration by parts. It depends on a single time variable ti. Using two functional relations, (2.32) and 



dx 



— [Qi,iuix) - Q-i,iuix)] 

ITT 



( see Appendix 2 ) its Fourier transform can be presented in the form 



which leads to the Poisson-type integral, 
with the instantaneous kernel, 



krQ-lAuik±Ti) 
vQl,iu{k^Ti) 



dr2drj2K.„i{Ti;fi - ^2, 771 - ?72)p(ti, ^2, '72) , 



dvdk e"'-'i+''^^ 
(27r)3 ik\ 



krQ-l^iu{k±_T) 
vQlAv{kx_T) 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



The potential given by Eq. ( 2.44 ) coincides with the potential V^^'^^ of Eq. ( 2.18 ) of the r-static source. 

Therefore, this term represents the instantaneous distribution of the potential at the moment ri , corresponding to the 
charge density taken at the same moment. Next, we have to recall that the charge density p{k, v, t\) in Eq. ( 2.44 ) still 
includes an arbitrary constant p^ik^v), which may be interpreted as the charge density at r = 0. This constant has 
appeared because we used only the conservation ( 2.24| ) of the constraint (which is the consequence of the equations 
of motion) and did not used the Gauss law explicitly. Now we can see that imposing the constraint indeed affects 
only the static potential of the charge distribution and puts it in agreement with Gauss law. If the initial data allow 
one to put po = p{t — 0) = 0, then it immediately solves two problems. First, the conservation of the constraint 
just duplicates the Gauss law, and the latter can be used to remove the unphysical degrees of freedom without any 



reservations. Eventually, it allows to fix the gauge AJ = ^ completely (see Sec. [II A). Second, it becomes possible 
to eliminate the charge density p completely, and to replace it by the components of the current. Replacement 
follows an evident prescription, 

dp 

P{Tl,r]2,r2) ^ dT2—=-l T2dT2[ksj'' {t2, 7^2, r2) + iyf{T2, 1^2, f2)] , 

Jo 0T2 Jn 

and leads to the standard form of the A'^™'**) representation (an artificial contribution of any po would correspond to 
the recognizable static pattern in the longitudinal part of the propagator and is easily handled), 

T2dT2 I dri2d'^r2^iZ''*\Ti;T]i - 772, n - r2)j'"(T2, 772, ^2) , (2.47) 



with the kernel given by the formula. 



A(inst) I 



dvd^k 
(27r)3fc2 



krQ-lauik±Ti) 
vQl.iv{k^Ti) 





■fc/ 


I 


V 



„w(?)i-?72)+ife(n-f2) 



(2.48) 



Two parts of the propagator, given by equations (2.42) and (2.48) can be combined in one elegant formula for the 
propagator of the field ) = A^-^^ + A^™"*) , 



A 



(long) ( _, N 

im (n;?7i -f]2,?-i -r2) = 

/■ dvd^k krks{Q-\^iAkx_T2) - Q-\,iu{k\_Ti)\ krV{Qx,iv[kx_T2) - Q-\.iv(k^T{)\ 

J (27r)3A;5_ [ '^ks[Q-i,iu{k±T2) - Qi,iu{k±Ti)] i^^[Qi,iu{k±T2) - Qi,iu{k±Ti)] 



(2.49) 



Im 
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This expression will be used for practical calculation of quark self-energy in paper [IV] , where it will be transformed 
into the mixed representation. The coordinate form of A'^^ that reveals its causal properties, is analyzed in Sec. |^ 
Eqs. ( 2.39| )-(2.42) and (2.47), (2.48) present the propagator in a split form. Different constituents of this form 
are as a preliminary identified as the transverse, the longitudinal and the instantaneous parts of the propagator. It 
would be useful to learn if the same kind of splitting is possible for the current itself. An affirmative answer (as in 
the cases of the Coulomb and radiation gauges) would be helpful for the design of the perturbation theory. To answer 
this question, one should substitute the different pieces of the solution into the left hand side of the original system 
of differential equations. This leads to the following expressions for the Fourier components of the three currents: 



'fcis-l,i,.(fc^T) 



dp 



(2.50) 



-( 




m 




V 





[Q-lMk^T){k^f + kyf) + Ql,zuik^T)f' 



(2.51) 



VT^^Ql,iAk^T) 





1 


dr) 





krSl,iu{kl_T) 

vk\s^i^iy{k^T) 



dp 
97 



(2.52) 



Provided that the current is conserved, these three currents, added together, give the full current on the right hand 
side of the system. Therefore, the solution is correct. However, none of these three currents carries any signature 
of being longitudinal or transversal in the usual sense. None of them has zero divergence since the operator of the 
divergence does not commute with the differential operator of the system. No desired simplification is possible in our 
case. 

In fact, the above splitting of the potential has no real physical meaning. To see it explicitly, let us find the 
divergence of the electric field, divE = dm£"^ [again, for brevity, in the Fourier representation]: 

,2 

(2.53) 



divE(*'^)(fc, v-t) = i[Q-i,Uki^T) - QiMk±r)]{iyT^f - ^{k^f + fc^j^)) , 



divE(^)(fc, ly; r) 



+ kyjy)Q-i^,^{k±T) + vfQl,^^,{k^T)] , 



(2.54) 



divE(™^*)(^, v-t) = p{k, v-t) - t{T^Q-iMk±r) - ^Q,^,^{k^T))[{k,f + kyf) + vf] . 



(2.55) 



Only the divergence of the true retarded component of the field E*^**") turns out to be zero. The term which prevents 
the divE'^*''^ from being zero is due to the non-symmetry of the propagator. A'''' ^ A'"''. It appears when the 6*- function 
in Eq. ( 2.39D is differentiated with respect to Hamiltonian time r. This term is vital for obtaining the expression that 
obeys the Gauss constraint, divE(fc, v\ r) — p(k, i/; r). 

The known examples, when the transverse and the longitudinal fields are separated at the level of the equations 
of motion, are related to a narrow class of homogeneous gauges. The impossibility of a universal separation of the 
transverse and longitudinal currents thus appears to be a rule rather than exception. It reflects a general principle; the 
radiation field created at some time interval has the preceding and the subsequent configurations of the longitudinal 
field as the boundary condition. The dynamics of the longitudinal field falls out of any scattering problem in its 
^-matrix formulation. However, this dynamics is, in fact, a subject of the QCD evolution in the inelastic high-energy 
processes. 
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D. Initial data problem in the gauge = 



We obtained the expression for the (retarded) propagator as the response function between the "external" current 
and the potential of the gauge field. We must also verify that the same propagator solves the Cauchy problem for 
the gauge field. This can be easily done by presenting the initial data at the surface t = tq in the form of the source 
density at the hyper-surface t = tq, i.e., 



lr{r2) = V^iRg""(ro)[5'(r2 - ro)A„,{r,i^) + 5{t^ - To)A'^{r,n)] , 



(2.56) 



where A„i (r , 77) and A'^ (r , 77) arc the initial data for the potential and its normal derivative on the hyper-surface 
r = To. Usually, it is assumed that the external current vanishes for r < To. Substituting this source into the 



Eqs. (2.39), (2.41) and (2.45), and taking the limit of ti — ^ To, we may verify that the standard prescription for the 
solution of the initial data problem. 



d 

(fr2 dri2 I\im{xi,X2)^A"\x2) 
OT2 



(2.57) 



holds with the same propagator /S.im{xi,X2) that was used to solve the emission problem. For example, in the limit 
of T — > To, the 77-component of the vector potential is a sum of three terms, 



4*'-)(to + 0) = ^ { [Rfl{k^To)H\l\k^To) - R^^l{k^ro)H\l\k^To)][uA^ - klA,] 
^Tov[R^^l{k^To)R%,{k^To) - R^'Uk^To)R%,ik^To)]A'^ } , 



(2.58) 



— — lyk — ly — 

9i,ii^ik±To)A^ + ToQ-i.ity{kj_To)A'. is_i,,v(fc_LTo)A,, H Qi,i„{k^To)A' ] , (2.59) 

To To 



A {inst 
V 



{to + 0) = ■^QiMk±To)[ToA'^ + —A'^] 

To 



(2.60) 



where we have denoted: A^ = k^A^ + ^y-^y. Eq. (2.6C) follows from Eq. ( ^.45| ) and takes care of the consistency 
between the charge density at the moment tq and the initial data for the gauge field. Using relations ( 2.32| ) and (2.44) 
and adding up Eqs. (2.5S)-( p.60D we come to a desired identity, A,,(to + 0) = A,,. 

When the initial data Am {r, rj) and A'^ (r, 77) correspond to the free radiation field, then only the part of the full 



propagator, A\^\xi,X2), "works" here, and only Eq. (2.58) may be retained. The other two equations acquire the 



status of being constraints imposed on the initial data. Since the current is absent, we have A*^-^^ = on the left hand 
side of the Eqs. (^!59|) . Then the right hand side confirms that the kernel JC is orthogonal to the free radiation field 
modes. Since the charge density p vanishes, we have A^™**^ = 0, which is equivalent to Gauss' law for the free gauge 
field. The two transverse modes already obey these constraints. This fact provides a reliable footing for the canonical 
quantization of the free field in the gauge A'^ = 0. Indeed, the Riemann function coincides with the commutator of 
the free gauge field. It can be found via its bilinear decomposition over the physical modes. Thus, one can avoid 



technical problems of inverting the constraint equations (see Sec. III.) The longitudinal part of the propagator will 



be studied in details in Sec. IV 
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E. Gluon vertices in the gauge = 



The terms proportional to the first and the second powers of the couphng constant in the classical wave equations 
may be viewed as the external current and allow one to define the explicit form of the 3- and 4-gluon vertices. One 
should start from the solution of the Maxwell equations, 

= J d^xAt:iiz,,x)y/^^J^ix) , (2.61) 

with the color current of the form 

-9^V^)fabcfc<ihe'''{x)g'^\x)A\{x)At{x)Al{x) . (2.62) 

In perturbation calculations, every field A{x) in the RHS of this expression is a part of some correlator A(x, zjv). 
The components of the metric depend only on the time r while the derivatives affect only the spatial directions 
m" = (r, 77). Moreover, in these directions, all the gluon correlators depend only on the differences of the coordinates 
and can be rewritten in terms of their spatial Fourier components. After symmetrization over the outer arguments 
zn, one immediately obtains, 

Kl,c"bl,P2,P3;r) = -ItUc Sip, +P2 +P3) [g'"(p2 -P3)' +g"'(p3 -Pi)' +g''(pi -P2)T , (2.63) 

where p" = g^^pk, and the components of the momentum in the curvilinear coordinates are equal to pk — {Px,Py, i^)- 
The four-gluon vertex has no derivatives and is the same as usually. 



III. QUANTIZATION 



The second quantization of the field has several practical goals. We would like to have an expansion of the operator 
of the free gluon field like 

A{x)= /d2fcd^[cA(i^,fc)Fj^^^/a;)+4(^,fc)V;(^^);(a;)] , (3.1) 

with the creation and annihilation operators which obey the commutation relations 

[cx{iy,k),cl,{iy',k')]^6xyS{iy~iy')S{k-k'), [cx{'y,k),cy{iy' ,k')] = [c[{iy,k),c{,{iy' ,k')] = . (3.2) 



Once obtained, the commutation relations (3.2) allow one to find various correlators of the free gluon field as the 
averages of the binary operator products over the state of the perturbative vacuum and express them via the solutions 
V^-' (x) and F*--' (x). For example, the Wightman functions, 

iAio,.,(x,y) = {0\Aix)A,iy)\0) = E / = ^^oiMv,^) , (3-3) 

serve as the projectors onto the space of the on-mass-shell gluons and should be known explicitly in order to have 
a good definition of the production rate of gluons in the final states. The Fock creation and annihilation operators 
are also needed in order to define the occupation numbers and to introduce the gluon distributions into various field 
correlators defined as the averages over an ensemble. With these two Wightman functions at hand, one immediately 
obtains the expression for the commutator of the free field operators, 
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Ao,^J (x, y) = -i{0\ [A,{x), Aj (y)] |0) = Aio,^ (x, y) - Aoi,,^ (x, y) , (3.4) 

which should coincide with the Riemann function of the homogeneous field equations. The program of second quan- 
tization does not reveal any technical problems if we give preference to the holomorphic quantization which is based 



on commutation relations (3.2) for the Fock operators. However, if we prefer to start with the canonical commutation 
relations for the field coordinates and momenta, then one should postulate them and derive ( |3.2| ) as the consequence. 

The way to obtain the canonical commutation relations in cases of the scalar and the spinor fields is quite straight- 
forward. For the vector gauge field, we meet a well-known problem, viz., an excess of the number of the components 
of the vector field over the number of the physical degrees of freedom. For example, in the so-called radiation gauge, 
A° = and divA = 0, we write the canonical commutation relations in the following form W, 

[^.(x,t),i?,(y,t)] =<5-(x-y)= j ^ (^J,,- - ^) e^^^^^) , 

[A,{^,t),A,{y,t)\ = [S,(x,t),£;j(y,t)] = , (3.5) 

thus eliminating the longitudinally polarized photons from the dynamical degrees of freedom. The function plays a 
role as the unit operator in the space of the physical states. Here, z, j = 1, 2, 3 and the number of relations postulated 
by equations (|3.5|) apparently exceeds the actual number required by the count of the independent degrees of freedom. 



A = 1,2, of the free gauge field. The Fourier transform of the function 5*J is easily guessed because the basis of 
the plane- wave solutions is very simple [Q , and it can be obtained rigorously by solving the system of the constraint 
equations A similar guess or procedure in our case is not so obvious. We have the gauge condition A'^ = as 

the primary constraint and Gauss' law as the secondary one. The latter can be resolved in a way which allows one 
to exclude the 77-components of the potential and the electric field from the set of independent canonical variables. 
Thus, only x- and y-components are subject to the canonical commutation relations. To resolve the constraints, one 
needs the integral operators with the kernels built from the solutions of the Maxwell equations in the gauge A'^ — Q. 



Therefore we shall proceed in two steps. In section [HB, we shall sketch the results for the Wightman functions 



(3.3). These, will be used for the explicit calculation of the free field commutator (3.4) in section IHC and for the 



study of its causal behavior. 

A. Fixing of the gauge A"^ = 

Only the independent components of transverse fields are the subject for quantization. In order to eliminate the 



extraneous degrees of freedom we have to incorporate the Gauss law. In Sec. HC, the latter was shown to be a 
consequence of the equations of motion only when the evolution begins with the zero charge density at r = 0. Are 
these initial data physical or do they mean that the QCD evolution begins from nothing? Addressing the hadron 
collisions, the question can be asked more specifically: do we really need any resolved ad hoc color charges (dipoles, 
quadrupoles, etc.) to initiate the color interaction? The answer is negative for three reasons. First, as it has been 
demonstrated in paper [II], at t — + the states of wedge dynamics are widely spread along the null-planes with zero 
density at t = 0. Second, the nucleus which fluctuate into the expanded state of wedge dynamics is colorless even 
locally. Third, as it will be demonstrated in paper [IV], the interactions which are the strongest at very early times, 
are magneto-static by their nature. The collision of two nuclei is more likely to begin with the magnetic interaction 
of the color currents in the locally color-neutral system than with the electric interaction of color charges. 
In Abelian case considered here the gauge transformation is 



13 



(3.6) 



Since we have At-{x) = 0, we must also have x{^) — x{'^t,fi)- The boundary condition, A^(r = 0,77) = 0, cannot 
be altered by the gauge transform (3.6). Therefore, we must require that x(^) = x(^)- Hence, the allowed gauge 
transform is reduced to x- and y-components of the vector potential. 



The Gauss constraint (2.12), 



A'^{x) = Ar{x) + drxin) ■ 



d^drA^{T,T]) = -T'^drAr{T,'n) - T^f{T,r]) 



(3.7) 



is a hyperbolic differential equation for the function Ay^ (r, 77) which can be integrated (the Riemann function of this 
equation, R{T,ri) — 1). With the boundary conditions, ^^(0,77) = and A^{t,— 00) — 0, this equation has a unique 
solution. 



^r,(r, 77) 



d77 / dT[~T drAr(j,ri) - T fir, if)] . 



(3.8) 



The residual gauge transform (3/7) changes only the integrand of Eq. (3.8) 

drK{r) ^ drArir) + ^^x{rt) , 

As a consequence of the boundary conditions, the transverse divergence of the field must vanish at r = 0, 9^ ^^(0, 77) = 
0. Therefore, we must also have 

A±x(rl) = , 

x(x,y) must be a harmonic function. Demanding, that xi^jU) vanishes at |rj| — > 00, we find that x(a:,y) — 0. The 
gauge = is fixed completely. The Gauss law can be unambiguously used to eliminate A,j from the list of canonical 
variables. 



B. Gluon correlators in the gauge A"^ — 



In this section, we shall write down components of the field correlator Aio,ij (a:, y) in the curvilinear coordinates 
u — {T,ri,f). We shall denote their covariant components as Aio,i/c(Mi7 1*2). Later we shall transform them to the 
standard Minkowski coordinates and find the correlators of the temporal axial and the null-plane gauges as their 
limits in the central rapidity region and in the vicinity of the null-planes, respectively. The most convenient (for this 
purpose) basis consists of the transverse modes v'^^K The mode v^^^ gives the following contribution to the correlator 
Aio,ife, 

iA^l) n 9^- f ( ^zhl\piHri-r2)„-ik±Ticosh{B-m)+ik^T2Cosh(e-ri2) (OQ\ 

»^io,..(i-^J-y_^yy (^(^d.. —jee . (3.9j 

Realizing that de/2 = dk^/2k°, we recognize a standard representation of this part of the correlator in terms of the 
on-mass-shell plane waves decomposition. 

The second part of the correlator is determined by the mode v^^^ and has the following components: 

A(^„),,(l,2) = ^1 _g_^e'Mn-.^.);^(0,,^,^^)/*(0,,^,^^) , (3.10) 
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Aro:,,(l,2) = -* / - 



(27r)3 



(3.11) 



/oo in r rj'^h p'^k{f\—r2) 
^YJ (2^)3 fc2 f2i0,r„r„)f;{e,r„7j,) . 

One may easily see that all components of Aio(l, 2) vanish when either ri or T2 go to zero. 



(3.12) 



C. Causal properties of the field commutators in the gauge = 

Causal properties of the radiation field commutator may be studied starting from the representation (pM. Using 



Eqs. (3.9) and (3.10) we may conveniently write the contribution of the two transverse modes in the following form, 



(27r)3 



(3.13) 



zA(^),(1,2) = -, 



d k kip k 



(27r)3 k\ 



d9 



cosh 2t] 



sinh^ 9 + cosh^ rj 



(sinfcj^<i>i2 - sin/cj_$i + sinfcj^(f>2) , (3.14) 



where we have introduced the following notation: 2ri = rji — ri2 , r ^ — r2 , $i = ti cosh(6' — rji) , $12 = $1 — $2 
The sum of (|A3|) and (|I|) can be rearranged as follows, 



iAo,^.(l,2) = 



(fk de 
72^^ 



ikr 



k k 

6rs sinfc_L$i2 + -^-i2--[suik_\_'^i — sinA:_L$2] 



+fc,.fcsCOsh(77i - 772) / dr' j dr" sin[A:±r' cosh(6l - yy) - fc^r" cosh(6l + 77)] . (3.15) 

JQ Jo J 

Rewriting the integration d^k dO into the three dimensional integration d^if^li^l jj^ Cartesian coordinates, the first 

d'^k dO . . ^ sign(ti - 



integral in ( p. 151) 



Z?o(l,2) = 



-e""" sinfc^$i2 = 



27r 



(3.16) 



is easy to calculate and to recognize as the commutator of the massless scalar field. It differs from zero only if the line 
between the points xi and X2 has a light-like direction. We integrate the first and the third terms in the integrand 
of the Eq. ( plsj ) in this way. To reduce the two integrals in the second term to the same type, we must exclude the 
factor 1/fc^ using the fundamental solution of the two-dimensional Laplace operator. 



krks 
kl 



drds I ^ln|e-r1e'^^« 



(3.17) 



After that, we arrive at the final result, 

Ao,.,(l, 2) = ~5rsDo{l, 2) - cosh(r;i - r;2)a,.9. / ' / ' dr2i?o(l, 2) 

JQ JO 



+ drds 



(27r) 



(3.18) 



From this form, it immediately follows that the commutator of the potentials vanishes at ri = r2. An even stronger 
result is found for the commutator of the two electric fields, 
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52 



dTldT2 



iAo,r.(l,2) = 



-Sr 



92 



9ti(3t2 



cosh(77i — 772) 



92 



(3.19) 



This commutator vanishes everywhere except on the light cone, in full compUance with the microcausahty principle 
for the electric field which is an observable. However, this does not happen for the commutator of the potentials since 
they are defined non-locally. It vanishes neither at space- like nor at time-like separation because the line of integration 
which recovers the potential at the point X2, in general, intersects {e.g. at some point xs) with the light cone which 
has its vertex at the point xi, and the commutator of the electric fields at the points xi and is not zero. 

Similar results take place for the commutator of the 77-components of the potential and the electric field. The field 
commutator, 



[E^{l),E,,{2)] = 



92 



9ti9t2 

is entirely causal, while the commutator of the potentials, 



zAo,r,^(l,2) = -iV\Da{l,2) 



[^^(1), ^^(2)] = iAo,^^(l, 2) = [ ' ndn / ' T2dT2Do{l, 2) 

Jo "'0 



(3.20) 



(3.21) 



does not vanishes at space-like distances, except for ti ~ T2. Finally, the formally designed commutator between the 
r- and 77-components of the electric field (the two observables). 



[Er{l),E^{2)] 



92 



-iAo^r„(l,2) 



92 



-^0(1,2) 



(3.22) 



9ti9t2 ' ' dx^dri\Ti 
is entirely confined to the light cone, while the commutator of the potentials (which are not the observables), 

[A,(l),A^(2)]=zAo,,^(l,2) = -y^ ^2d^2__Z?o(l,2) , (3.23) 

does not vanish at the space-like distance, even at ti — T2. This result, however, is not a subject for any concern since 
the potentials are defined non-locally and commutation relations for the electric and magnetic ( cf. (2.14) ) fields 
are reproduced correctly. Moreover, we have argued above that the //-components of A and E are not the canonical 
variables since the constraints express them via the x- and y-components. 

The "acausal" behavior of the Riemann function, Aq''(1,2), may cause doubts whether the gauge A'^ — allows 
for meaningful retarded and advanced Green functions which, by causality, should vanish at space-like distances. For- 
tunately, this anomalous behavior appears only for the gauge-variant potential; the response functions for observable 
electric and magnetic fields are causal. This can be easily seen, e.g., from Eqs. (2.21), (2.25) and ( p. 26 ), which are 
the usual inhomogeneous relativistic wave equations for the various physical components of the field strengths £ and 
B. 



D. Canonical commutation relations in the gauge A'^ = 



A proof of the commutation relations ( |3.2D for the Fock operators follows the standard guidelines [[7|. First, the 
creation and annihilation operators are defined via the relations, 

c;,(7.,fc) = ivf,A) = tg^^ / d3x[FW*(x)A,(x,r) - T/ . ,(:r)A,(x, r)] , 



cliiy, k) = {A, V^;^^) = 7g*-'' / d3x[A,(x, T)v%i^) - A.(x, r)V^^^^{x)] . (3.24) 
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This results in the following expression for the commutator, 



.(A)* .(A') 



+ [M^,T),Ai{y,T)] V,i,,{x)V,r^,.Jy) + [M^,T),Ai{y.T)] V^^y{x)vl^]^{v)] . (3.25) 

Most of the terms in the second hne vanish due to the commutation relations. Next, we rely on the following guess 
about the form of the commutator, 



(3.26) 



A=1.2 



which leads to the proper equal-time commutation relations for the independent canonical variables. Finally, explicitly 
using the orthogonality relations for the eigen- modes we immediately obtain the commutation relations (|3.2D. 



IV. LONGITUDINAL PROPAGATOR AND STATIC FIELDS 



In this section, we shall find the explicit expressions for the kernels (2.42) and (2.46) which represent the longitudinal 
and instantaneous components of the gauge field produced by the "external" current j'^. The calculations are lengthy 
and their details can be found in Appendix 3. Here, we present only the final answers. 

The components of the longitudinal propagator a[^^ (t2 , r, rf) are already obtained in the form of the three- 



dimensional integrals (2.42). A^^^ depends on the differences of the curvilinear spatial coordinates, r ~ ri — r 2 
and ?7 = r]i — r]2, but not on the difference of the temporal arguments ti (of the field) and T2 (of the source). 
Introducing the shorthand notation for the distance in the (a;y)-plane, r±_ |?^, and for the full distance 
i?2 = R{t2) ~ [{ri — ^2)^ + T2 sinh^(77i — 772)]^^^ between the two points of the surface T2 = const, we obtain: 



(t2 - r±_) 
An 



1 - 



T2 cosh 7] 



R2 



6r 



2x^x'^ T2 cosh 77 



^(L) ^ 0{T2-rj_) x' T2sinh ?? r| 



Ai^) 



y(5(rO(5(r;) 



0{r2 - r^) 
47r 



i?2 

^ 7] coth — 1 
sinh^ 77 



9{t2 — r±) x'"" r| sinh77 



T2 cosh 77 / 
i?2 sinh^ 77 \ R2 



4tt 7-^ 

2 cosh 77 



sinh 1 7; I 



(4.1) 



where L2 — L{t2) ~ ln[(T2 sinh I77I + R2)/r\}^. By examination of Eq. ( ^.4l| ), one may see that after the replacement 
of T2 by Tl, the same kernel, A|^-'(ri, r*, 77), determines the components E"™ '(n) of the longitudinal part of the electric 
field via the components j'"(ri) of the current at the same time. These propagators do not respect the light cone, but 
have a remarkable property that the longitudinal fields at the surface of the constant proper time r do not exist at the 
distance from their sources that exceed r. This establishes the upper limit for the possible dynamical correlations 
between the longitudinal fields in the (a;?/)-plane. 

One more representation of the longitudinal propagator is interesting at least in two respects. First, for the practical 
calculations in the paper [IV], we shall need the longitudinal part of the propagator in the mixed representation. 



^Im"^' (ti , T2 ; 7^1 — r]2,kt), which can be shown to be 



kt cosh 1 7; I 



^-tkt sinh |J7|^^ 



— tanh(Q; + 2 ) tanh(Q; — — ) [ sin ktTi — sin ktT2 ] 



(4.2) 



\[lony] _ _ll 



kt cosh 1 77 1 



-tkt sinh \ri\ ^2^^ 



da sin fc^Ti — sin ktT2 — kfTi cos fctTi + fcf 72 cos kiT2 



2ti 



hi cosh^(a + ^) cosh^(a - f ) 



(4.3) 
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where the last integral terms in Eqs. ( [l.2| ) and (4.3) provide that the longitudinal part of the field obey the boundary 
condition we imposed at r — )■ 0. These terms cancel with the similar terms in the radiation part of the retarded 
propagator A[^e(] (ti, r2). Eventually, the last fact guarantee that acausal terms (which were the subject for concern 
in the course of the canonical quantization in Sec. Ill D) do not contribute to the dispersion equation that we derive in 
the next paper. Second, the first very simple by its structure contact term in the component d|,'^"^' of the propagator, 

A[~l(r2, n; - m-,kt) - -^^-^^(ry) , (4.4) 

appears to be the only part of the enormously complicated full retarded propagator aJ^*' which significantly con- 
tributes the amplitude of the forward quark-quark scattering at the earliest stage of the collision. 
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APPENDIX A. Modes of the free gauge field 

Here, we shall obtain the complete set of the one-particle solutions to the homogeneous system of the Maxwell 
equations with the gauge — 0, that is, Eqs. ( ^.10 ) and (2.11). This gauge condition explicitly depends on the 



coordinates, thus introducing effective non-locality in the path integral that represents the action. Therefore, it 
becomes impossible to invert the differential operators using the standard symbolic methods. An explicit form of the 
one-particle solutions becomes necessary in order to find the Wightman functions of the free vector field, to establish 
the the explicit form of the field commutators, and to separate the propagators of the transverse and the longitudinal 
fields. It is natural to look for the solution in the form of the Fourier transform with respect to the spatial coordinates, 

A,{x) = / di^ f d'^k e"''^e'''^ A,{k,v,T) . (Al) 



Then, the system of second order ordinary differential equations for the Fourier transforms takes the following form: 

[Td^+dr + — +Tkl]A,^{k,l^,T)~Tk^kyAy{k,l',T) A^{k , 1^, t) = , (A2) 

— * jy"^ — » jyk — ' 
-Tk^kyA^{k,v,T) + [Td'^+dr^ + Tkl]Ay{k,iy,T) ^A„(fc, r) = , (A3) 

T T 

- '^A^ik, V, r) - ^A^(fc, V, t) + [^92 _ 1 5^ + -k\\A,{k; r) = . (A4) 
r r r T 

In this form, the system is manifestly symmetric and self-adjoint. An additional equation of the constraint reads as 

C{k,V,T)^^vdrAr,+Tdr{kxA^{k,V,T)+kyAy{k,V,T)\=^ . (A5) 
T 

Let us rewrite the homogeneous system of the Maxwell equations in terms of the variables 
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^ = d^A.^ + dyAy, = dyA^,~d^Ay, andA = A^. (A6) 

One immediately sees that the equation for the Fourier component ^'(fc, t) of the fongitudinal magnetic field ^'(a;) 
decouples, 

[dl + -dr + ^ + fci]*(fc, v,t)={) . (A7) 

T 

Then, the other two equations of motion take shape, i.e, 

[T^dl + Tdr + A'^kA^) - *^fciA,- ,(r) = , (A8) 

[92 _ + fci]A,, ,,(t) + iv^j:^^{t) = . (A9) 
The additional constraint equation can be conveniently rewritten as 

— ' lU 

C(k,V,T) = -drki^^{T)+Tdr^^^^^{T)=Q . (AlO) 

This is an independent equation. However, the conservation of the constraint along the Hamiltonian time t is 
a consequence of the equations of motion, and it can be employed to obtain the independent equations for the 
components of the vector field. This is easily done in terms of the auxiliary functions, 

nA^) = ^'^k,ui^) and = r-iA,;^(T), (All) 

which are directly connected to the "physical" components of the electric field, 5™ = \f—ZZ™'^Ai : 

1 



dr[dl + -dr+^-^ + ki]v^^,{r)^Q , (A12) 



dr[T''dl + Tdr + + fcir^Ja^- ,^ (t) = , (A13) 
As a result, we see that the functions (^^ ^(r) and a-^ ^{t) obey inhomogeneous Bessel equations, 

[dl + \dr+'^ + kl]^-^^^{T)^Pc^ , (A14) 

[^r + ^9, + ^ + fci]a,-^(r)=T-2c, , (A15) 

where and Ca are arbitrary constants. We may now cast the solution of these equations in the form of the sum of 
the partial solution of the inhomogeneous equation and a general solution of the homogeneous equation, 

^u,,{T)=aH^X{k^T)+a*H^Xik^r) , (A16) 
nA"-) = + c*H'^X{k^T) + c^s.^k^r) , (A17) 

a,, ^(r) = iH^ik^T) + i*HX{k^T) + c„s_i,,,(fc^r) , (A18) 
where s^^^{x) is the so-called Lommel function 1^,^. 
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Furthermore, it is useful to notice that the system of the Maxwell equations (2.10) and (2.11) also has an infinite 
set of the r-independent solutions of the form 



Wi{r],r) = d^xiv^r) 



(A19) 



where x is an arbitrary function of the spatial coordinates rj and r. Thus, they are the pure gauge solutions of the 
Abelian theory, compatible with the gauge condition. 

In order to find the coefficients one should integrate Eqs. ( A17 ) and (AlS) with respect to the Hamiltonian time 
T, thus finding the functions $ and A. Next, it is necessary to solve Eqs.(A6) for the Fourier components of the 
vector potential and to substitute them into the original system of Eqs. (A2)-(A5). Using functional relations from 
Appendix 2, one obtains that c+ i/j ^ and Ca — J^c^ = . 



One of the solutions, (already normalized according to Eq.(^.g|)) is found immediately: 



25/27rfci 







(A20) 



Initially, the components of the vector mode l^'-^-' appear in the following form required by the convergence of the 
integral, 

\ -i?fL,J/fcj.r|s) -zi/[e'^^/Vsinh(7ri^/2)] 



However, it can be gauge transformed to the more compact form, 



kr-vR 



(2) 



-R^l,,{k^T\s) 



(A21) 



The third solution (the last one by the count of the non-vanishing components of the vector potential in the gauge 
= 0) has the following form, 



krQ-l.iv{k^T) 
vQl,iv{k±_T) 



(A22) 



The modes V^^'^ and F'^) are the normalized solutions of the Maxwell equations. They are orthogonal and obey 
the normalization condition. 



(Fi'''\ yi'^')) = 6{v - v')5{k - k'), (T/i'\ yi'O = 

k.v k'v' k.u k.v 



T/(2)^ 



(A23) 



which can be easily verified by means of Eq. (B5). The norm of these solutions is given by the Eq. (2.9). A 
normalization coefficient of the mode V'^^^ is not defined, as this mode has a zero norm. It is also orthogonal to V^*^^^ 
and 



(yi^) yi^) . ^ . (1) ^(3). ^ ,ym (3). ^ 



k.u k^u 



k^u k,u ' 



(A24) 



Thus this mode drops out from the decomposition of the free gauge field. 

The conservation of the constraint can be obtained as a consequence of the Eqs. ( A12 ) and ( A13 ) in the form, 

^^AnA-") + ^H^-^y^ = ^drCj:u{T) = , (A25) 
which reassures us of the consistency between the dynamic equations and conservation of the Gauss' law constraint. 
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One can explicitly check that the modes V^^^ and T^^^^ obey the constraint equation (AlC), which expresses Gauss 



law. The mode V^^^ does not. This mode corresponds to the longitudinal field which cannot exist without the source. 

APPENDIX B. Mathematical miscellany 

This appendix contains a list of mathematical formulae for the functions which appear in various calculations in 
the body of the paper and Appendix 1. The components of the vector field are expressed via two types of integrals. 
The first of them was studied in Ref. |^ : 

Rlili^lS) = J x^^Hi=\x)dx = x[{p. + V - l)i/(^")(a;)S,,_i,._i(x) - hI'\{x)S ^A^)] , (Bl) 

where S^^i, stands for any of the two Lommel functions, s^^i, or S^^i^ H;^- [ Whenever we omit the indicator |S), the 
function R\i}v{x\s) is assumed.] The second type of integrals, 

Q^iA^) = / x^'dxs^^Ax) , (B2) 
Jo 

is a new one. The functions R\i}u{x\S) are introduced as indefinite integrals. The preliminary choice of the lower 
limit and, consequently, the choice of which of the functions, ,y or S'^^^^jis used is motivated by the requirement of 
convergence and regular behavior. One can easily prove that 

(2\ f-2\ -p- 7ri//2 



sinh(7rz^/2) 



r\%AAS) - r\%AAs) = -jj::^,. (b3) 



We often use the following relation between Lommel functions 

Si,,.(fcxT) = 1 - z/2s_i_,,(fc^T) . (B4) 



From the integral representations ( Bl ) and (|B2| ) , it is straightforward to derive the functional relations 



R%Ak^T) + ^R^^l{k^T) = -^^H\i\k^r) , (B5) 
-i,t^ik±T) - (3i,i,y(fc±r) = y-;r-si,i^{kj_T) = T—s^i^i,y{k±T) . (B6) 



OT OT 



The Wronskian of the Hankel and Lommel functions, 

W{s,,.(x),h\^{x)} = -\r^I{x) , (B7) 



must be obtained in order to prove orthogonality of y'^' and V'^^^ To prove (B7), one should use the following 
representation for the Lommel function, 

s,,M = - Hl^\x)R[%ix)] , (B8) 



which follows from Eq. (Bl) and, consequently, 



s'lM^) = l^[H^l^'{x)R^^l{x) ~ H^'J (x)R^^l{x)] . (B9) 
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In order to prove relation ( 2.32| ) one should use the representation ( |B1| ) for the functions and the Wronskian 
of the two independent Hankel functions. The proof of relation ( 2.43 ) begins with replacing the functions r''_1\ by 
^^i\u means of Eq. ( B5 ) . The final result follows from Eq. ( B9 ) and (B6). 

APPENDIX C. Calculation of the longitudinal part of the propagator 

The kernels ( 2.44 ) and ( ^.46 ) of the longitudinal and instantaneous parts of the propagator are given in the form 
of the three-dimensional Fourier integrals dv(Pk. Here, we describe the major steps of the calculations which lead to 
Eqs. (13). 

We permanently use the following integral representation for the Hankel functions, 



^ J —OO 



(CI) 



which allows one to calculate many integrals by changing the order of integration. The Lommel function Si^ii, has a 
similar representation. 



SiAv{x) — X cosh u COS z^u e 



-X sinh u 



du 



Integrating by parts, and using Eq. (B4), we find the integral representation for S^i,i^, 



sai{vu) e 



■X sinh u 



du 



(C2) 



(C3) 



We start with the integral representation (B2) of the functions Q±i^ii, and perform an integration over ly. To compute 
the integrals from the function si^i, it can be conveniently decomposed in the following way. 



which allows one to find 



-7ri//2 



sinh(7rz^/2) 



(C4) 



Si.,^{k^T)e"''^dv = Trfcj^rcoshr/e"'^^"'"'^!''! 

) 

/OC 
vS-l^^^{k^T)e^'''^dv = iTrsignr^g-'^^^"'"'^ I"! 
-00 



(C5) 



The similar Fourier integrals from the function are calculated using the representation ( |Cl| ) for the Hankel functions 
and the integral. 



sinh(7ri^/2) 



dO- 



cosh 



This yields, for example, 



dve"''^h,^{kx_T) = 



2 — sin[fc^T cosh(0 — i])] 

00 cosh 6 



After integration over z/ we obtain the following integral for Ai-s , 

'■^^ dr 



d'^k krks ^^f:p 



(27r)3 kl 







Trfc^rcoshrye-'^^^'^'"'^!''! 



-00 cosh 



■ sin[fc^T cosh(0 ~ I i 



(C6) 



(C7) 



(C8) 



and similar integrals for the other components. The first term in this formula is calculated in the following way. After 
integration over r we continue: 
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coth \r,\ r ^Jo{k^r^)[l ~e~^^^' '^'"'^ I"!] = - Mi coth hi 
Jo k± Ait 



T2 sinh \rj\ + y T2 sinh 77 



(C9) 



To work out the second term, one should introduce = k± sinh 6* and fco = k± cosh 6* = |k| and change (Pkd9 for the 
three-dimensional integration d^h. With t = r cosh 77, r ~ {x, y, r sinh 77), this leads to 

T cosh 77 



„ drds n dr f d^k • , ^ drd. 



r'j^ + sinh^ 77 



An 



{rj_ — T2) coth I77I In 



T2 sinh I77I + \/rj_+T2 sinh 77 



+ 9{t2 -r^) In 



T2 



(CIO) 



Adding (|C9| ) and ( |C1Q| ), we obtain the first of the equations (4.1). 

APPENDIX D. Gluon correlators in the central rapidity region and near the light wedge 

In this section, we compare the correlators of the gauge = with the similar correlators in the three other 
gauges, A" = 0, — and — 0. We shall start with the simplest on-mass-shell Wightman function A^'q . These 
type of correlators, Aq{', Aq" and A^'' share the same polarization sum of the free gauge field. They correspond to 
the densities of the final states of the radiation field and are important for various calculations. The same polarization 
sum appears in expressions for the transverse part of the propagators, Ajf^j, A[^^^, Aqq and A^'^. For our immediate 
purpose we shall include the projector d'"' of the gauge = to the formal Fourier representation, 

d^k 



iD'^^ {xi,X2) 



rd''"(fc;a;i,X2)e-*'=("i-"^\ 



(Dl) 



(27r)32fcO 

with the "extraneous" dependence of the Fourier transform on the time and spatial coordinates. This dependence 
disappear in some important limits. Therefore, we discover the domains where the wedge dynamic simplifies and 
describes the processes which are approximately homogeneous in space and time. These domains are: (i) the central 



rapidity region, 771^2 ^ 1 (or 2 ^ 0)i where the projector in the integrand of Eq. (Dl) is 



rf"'(fc,u) = 



(D2) 



ku (ku)^ ' 

with the gauge-fixing vector, = (1, 0, 0, 0), which approximately coincide with the local normal to the hypersurface 
T = const; and (ii), the vicinities of two null- planes, 77 — > ±00 (or x^ — > 0) where 



d^^-{n±,k) = -g^- 
with the null-plane vectors rij^ = (1, 0, 0, =Fl) . 



{kn±) 



Eqs. (3.9)-(3.12) almost fit our needs. In all three cases {x^ 0, k'^x'^ 3> 1, as well as 



(D3) 



0, krx+ > 1 and 



0, k^x 3> 1) the functions /i and /2 can be approximated by the following expressions: 



fi~i tanh(6' — ri)e 



-ik± T cosh(9—ri)-\-ikr 



k^'^ x^ k"^ x^ 



-ikx 



k^x 



k x"* 



-ikx 



k x+ 



(D4) 



/2 « ik±T 



-ik±T cosh{6 — ri)-\-ikf 

cosh(6' — 77) 



■ I 2 2 

= ik , T 



-ikx 



^ k°x" - k^x^ 



= 2ik^, 



-ikx 



k^x + k x+ 



(D5) 
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(Wc have omitted the time independent terms in /i and /2 which set the potentials of the mode v^"^^ to zero at r = 0. 
These kind of terms would correspond to the residual gauge symmetry and is not kept in the axial and the null-plane 
gauges as well. Thus, we cannot really claim the correspondence of the longitudinal fields between the wedge dynamics 
and these three dynamics.) 

Transformation of the correlator A''"(l, 2) to the Minkowski coordinates is carried out according to the formula, 

D^''ixi,X2) = <(a;i)g''(ari)A;„(wi,«2)g"'=(ar2)aj;(ar2) , (D6) 
where the matrix of the transformation is defined in the standard way, 

<(^)=a^' <(^)=^'> <(^)=^°. < = S: . (D7) 
These are the only components of the tensor af (a;) which participate in the transformation. In this way, we obtain, 

2) = a;?a;3A'"'(l, 2); D^^l, 2) = a;?a;°A'"'(l, 2) 
1)30(1, 2) = x^x^A""!!, 2); D^'{1, 2) = .t?xOA'"'(1, 2), 
^^'^(l, 2) = a;?A'"^(l, 2); £''■0(1, 2) ^ x3A"'(1, 2); 
2) = a;;A'"-(l, 2); £'''3(1, 2) = x^A'-''(l, 2), 

D'''{1,2) = . (D8) 

Every additional factor g^'^ = finds a counterpart which prevents singular behavior at r = 0. In the above 
approximation, the expression for the A''''(a;i, 0:2) component of the correlator has the form, 

J (27r)32fcO (fcOa;0 _ k^xDik^x^^ - k^xl) J (27r)32fc0 {k+x^ + k-x+){k+x^ + fc-^) ' ^ ' 

Therefore, in the limit of a;f 2 ~* we obtain that D^^,D^^ 0, while (P^{k,u) k^/k^, thus reproducing the 
corresponding components of the gauge = 0. The other components are reproduced one by one as well, and one 
can expect a smooth transition between the gauge of the the wedge dynamic and the local temporal axial gauge of 
the reference frame co-moving with the dense quark-gluon matter created in the collision. 
In the limits of a;f 2 ^ we obtain that 

1)00,1)03,1)30,1)33 ^ ^ ^ ifa;-^0, (DIO) 



Al- 



and 



Z)00, 2)03, ^30^ ^33 _ = ^ , if a;+ ^ . (DU) 

These limits, after they are found for all components, lead to the well known expressions of the projectors d'^"{n, k) 
in the null-plane gauge (the gauge = in the vicinity of a;+ = 0, and A~ = in the vicinity of x~ =0). Therefore 
we obtained an expected result; in the limit of the light-cone propagation, the gauge = recovers the null-plane 

gauges = and A~ =Q. 

Some remarks are in order. First, the concept of the structure functions relies heavily on the null-plane dynamics 
which essentially uses these gauges. For two hadrons (or two nuclei) we have two different null-plane dynamics which 
do not share the same Hilbert space of states. Now we have an important opportunity to describe both nuclei and 
the fields produced in their interaction within the same dynamic and the same Hilbert space. Second, one may trace 
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back the origin of the poles {ku)~^ in the polarization sums of axial gauges (uA) = and see that they appear in the 
course of the approximation of the less-singular factor [k± cosh(0 — rj)]^^ in various limits of the propagator of the 
gauge A'^ ^ 0. 

Further, contrary to the naive expectation that we obtain the gauge yl+ = at = and the gauge A~ = at 
x+ — 0, we obtained them in the opposite correspondence. First of all, let us notice that the result is mathematically 
consistent. Indeed, the gauge condition A'^ ~ may be rewritten in the form, 

^ ^{A'^e-'i + A-e") = . (D12) 

Thus the limit of ?/ oo {x^ — > 0) indeed leads to A^ — and the limit of ?/ — oo (a;+ 0) leads to yl+ = as 
the limiting gauge conditions. Recalling that 

A'^ = ^{A+e-''^ ~ A-e'i) A+e-" ^ -A-e'^ , (D13) 

we immediately realize that in the vicinities of both null-planes, the tangent component ~ 0. This fact has a very 
simple geometrical explanation; the normal and tangent vectors of the null plane are degenerate. Once = 0, we 
have = and consequently, = and A^ = at ?7 — ^ ±oo. This result naturally follows from the geometry of 
the system of the surfaces where we define the field states. These are subject to dynamical evolution in the direction 
which is normal to the hyper-surface. 



[1] A. Makhlin and E. Surdutovich, Phys.Rev. C 58 (1998) 389 (quoted as paper [I]). 

[21 A. Makhlin, Scenario for UHrarelativistic Nuclear Collisions: II. Geometry of quantum states at the earliest stage, hep 

ph/0007300| (quoted as paper [II]). 
[3] P.A.M. Dirac, Rev.Mod. Phys, 21, 392 (1949). 
[4] L.McLerran, R. Venugopalan, Phys.Rev. D49 2233, 3352 (1994). 
[5] H.Bateman, A. Erdelyi, Higher transcendental functions, McGraw-Hill, 1953 
[6] G.N. Watson, Theory of Bessel functions. 

[7] J.D. Bjorken and S.D. Drell, Relativistic quantum fields, McGraw-Hill. 
[8] S.Weinberg, The Quantum theory of fields, Cambridge Univ. Press, 1995 

[9] D.M. Gitman, I.V. Tyutin, Canonical quantization of fields with constrains. Springer- Verlag, 1990. 
[10] A. M akhlin and E. Su rdutovich, Scenario for Ultrarelativistic Nuclear Collisions: IV. Effective quark mass at the early 

stage. hep-ph/000730S (quoted as paper [IV]) 
[11] L.V. Keldysh, Sov. Phys. JETP 20 (1964) 1018; E.M. Lifshits, L.P. Pitaevsky, Physical kinetics, Pergamon Press, Oxford, 

1981. 



25 



